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s Abstract. We consider the effect of oblique EMIC waves on relativistic electron 

9 scattering in the outer radiation belt using simultaneous observations of plasma and 

10 wave parameters from CRRES. The main findings can be summarized as follows: 1. In 

11 comparison with field-aligned waves, intermediate and highly oblique distributions 

12 decrease the range of pitch-angles subject to diffusion, and reduce the local scattering 

13 rate by an order of magnitude at pitch-angles where the principle |n| = 1 resonances 

14 operate. Oblique waves allow the |n| > 1 resonances to operate, extending the range 

15 of local pitch-angle diffusion down to the loss cone, and increasing the diffusion at 
ie lower pitch-angles by orders of magnitude; 2. The local diffusion coefficients derived 
17 from CRRES data are qualitatively similar to the local results obtained for prescribed 
is plasma/wave parameters. Consequently, it is likely that the bounce-averaged diffusion 
is coefficients, if estimated from concurrent data, will exhibit the dependencies similar 

20 to those we found for model calculations; 3. In comparison with field-aligned waves, 

21 intermediate and highly oblique waves decrease the bounce-averaged scattering rate 

22 near the edge of the equatorial loss cone by orders of magnitude if the electron energy 

23 does not exceed a threshold (~ 2 — 5 MeV) depending on specified plasma and/or 

24 wave parameters; 4. For greater electron energies, oblique waves operating the \n\ > 1 

25 resonances axe more effective and provide the same bounce-averaged diffusion rate near 

26 the loss cone as field-aligned waves do. 
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2 ? 1. Introduction 

28 The flux of outer zone relativistic electrons (above 1 MeV) is extremely variable 

29 during geomagnetic storms. The competition between source and loss, both of which are 

3 0 enhanced during storm periods, determines the resulting relativistic electron flux level 

3 1 in the Earth’s outer radiation belt (RB) [e. g., Summers et al., 2004; Reeves et al., 2003; 

32 Green et al., 2004]. Usually, the flux falls by up to two or three orders of magnitude 

33 during main phase, and gradually increases over a period of a few days during storm 

34 recovery phase [e. g., Meredith et at., 2002]. Analyzing 256 geomagnetic storms during 

35 the period 1989-2000, Reeves et al. [2003] found that 53 % of storms lead to higher flux 

36 during the storm recovery phase in comparison to pre-storm levels; 28 % produce no 

37 change; and 19 % lead to net decrease in flux. The large electron flux decrease during 

38 the main storm phase is usually associated with a decrease of Dst when the relativistic 

39 electrons adiabatically respond to the stretching of the magnetic field lines caused by 

40 the formation of a partial ring current (RC) [Kim and Chan, 1997], and/or a drift 

4 1 out the magnetopause boundary [Li et al., 1997], and/or nonadiabatic scattering into 

42 the loss cone due to cyclotron interaction with electromagnetic ion cyclotron (EMIC) 

43 waves [Thorne and Kennel, 1971; Lyons and Thome, 1972; Summers and Thome, 

44 2003; Albert, 2003; Thorne et al., 2005] and/or whistler-mode chorus/hiss waves [e. g., 

45 Summers et al, 2007]. 

46 Precipitation of outer RB electrons due to resonant pitch-angle scattering by 

47 EMIC waves is considered to be one of the most important loss mechanisms. Recently, 
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48 data from balloon-borne X-ray instruments provided indirect but strong evidence for 

49 the ability of EMIC waves to cause precipitation of outer zone relativistic electrons in 
so the late afternoon-dusk MLT sector [Foat et at, 1998; Lorentzen et al., 2000; Millan 

51 et at, 2002], These observations stimulated theoretical and statistical studies which 

52 demonstrated that this mechanism of MeV electron pitch-angle diffusion can operate in 

53 the limit of strong diffusion, and can compete with adiabatic depletion caused by the 

54 D$t effect during the initial and main phases of storm [ Summers and Thome, 2003; 

55 Albert , 2003; Loto’aniu et at, 2006; Meredith et at, 2003]. 

56 Although the effectiveness of relativistic electron scattering depends strongly on 

57 EMIC wave spectral properties, unrealistic assumptions regarding the wave angular 

58 spread were made in previous theoretical studies. Namely, only strictly field-aligned or 

59 quasi field-aligned waves were considered as a driver for electron precipitation [e. g., 
so Summers and Thome, 2003; Albert, 2003; Loto’aniu et at, 2006]. The effect of oblique 
ei EMIC waves on relativistic electron scattering was recently discussed by Glauert and 

62 Horne [2005]. For prescribed plasma and wave parameters, considering the if + -mode 

63 EMIC waves, they calculated the equatorial diffusion coefficients and demonstrated that 

64 when a realistic angular spread of propagating waves is taken into account, electron 

65 diffusion at ~ 0.5 MeV is only slightly reduced compared with the assumption of 

66 field-aligned propagation, but at ~ 5 MeV, electron diffusion at pitch-angles near 90° 

67 is reduced by a factor of 5 and increased by several orders of magnitude at pitch-angles 

68 30° — 80°. As a result, EMIC waves should flatten the pitch-angle distribution. 

eg Thus, at energies of a few MeV, the assumption of field-aligned propagation 
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70 breaks down, significantly overestimating the pitch-angle diffusion coefficient at large 

71 pitch-angles, while underestimating the local diffusion rate at smaller pitch-angles by 

72 orders of magnitude. This is a very strong effect, so in contrast to [Glauert and Horne, 

73 2005 ], it is important to consider the impact of oblique EMIC waves on relativistic 

74 electron scattering using simultaneous observations for plasma/wave parameters, and 

75 to estimate the effect of bounce averaging. In the present study we calculate the 

76 pitch-angle diffusion coefficients using plasma and wave parameters observed by the 

77 Combined Release and Radiation Effects Satellite (CRRES) as reported by Loto’aniu et 

78 al. [ 2006 ]. 

79 This article is organized as follows: In section 2 we verify the pitch-angle diffusion 

80 coefficient calculations comparing our results with published results for both the 

si equatorial and bounce-averaged scattering rates. Then, using model wave spectra for 

82 H e + -mode EMIC waves with defined plasma parameters, we consider the effect of the 

83 wave normal angle distribution on relativistic electron scattering. In section 3 , using 

84 plasma/wave parameters observed by CRRES [Loto’aniu et al, 2006 ], we present the 
as results of our calculations and analysis of the local pitch-angle diffusion coefficients for 
as two selected wave packets. Finally, in section 4 we summarize the main findings of our 
87 study. 
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88 2. Equatorial and Bounce— Averaged Pitch— Angle Diffusion 

as Coefficients: Model Calculations 

An extensive statistical analysis of the EMIC events presented by Meredith et al. 
[2003] , showed that most of the cases when the minimum resonant electron energy fell 
below 2 MeV were associated with wave frequencies just below the He + gyrofrequency. 
So we take into account only the fJe + -mode EMIC waves in the present study. The 
model wave frequency spectrum is assumed to be Gaussian, 

5 2 (w)~expj — - — | > u lc < w < LOuCi (1) 

where, following Summers and Thorne [2003] and/or Albert [2003], ujlc — 
wuc — , u m — 3Qo+) — O.5f2o+, and flo+ is the gyrofrequency of 0 + . In 

our calculations, the wave normal angle distribution, g(9), is assumed to be a constant 
inside a specified region and zero otherwise. Below we consider the following three cases, 

Case A (field - aligned) : 0° < 6 < 30°, 150° < 9 < 180°, 

Case B (intermediate) : 30° < 9 < 60°, 120° <9 < 150°, (2) 

Case C (oblique) : 60° < 9 < 89°, 91° < B < 120°, 

where 9 is the wave normal angle. Note that the diffusion coefficient is a linear functional 
of the wave spectral density, and the sum of cases A, B, and C describe a situation 
when EMIC wave energy is evenly distributed over the entire wave normal angle region 
O°<0< 180° (we excluded the region near 90° because of Landau damping by thermal 
electrons [e. g., Thome and Home, 1992; Khazanov et al., 2007]). For benchmark 
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purposes, we calculate also the diffusion coefficients for a Gaussian distribution over 
x = tan 9 (0° < 0 < 15°) which has been used by Albert [2003], In each case, the wave 
amplitude is normalized to ensure 


J /wuc 

f d u ddB 2 (w) g (9) = 1 nT 2 . (3) 

uic Jo 

so Finally, to specify the ion content we follow Summers and Thome [2003], Albert [2003], 
si Meredith et al. [2003] , Loto ’aniu et al. [2006] , and prescribe the ion composition to be 
92 70% H + , 20% He + , and 10% 0 + (following [ Meredith et at, 2003] we call it a “storm 
93 time” ion composition). 

The results obtained using the relativistic version of the diffusion coefficient code 
of Khazanov et al. [2003] are shown in Figure 1. The first row shows the equatorial 
pitch-angle diffusion coefficients, the second row shows the corresponding resonance 
numbers averaged with the following weights: 


(n ( E,a )} 


n (g do; £ d9D n aa (u>, 8, E, a)) 


(4) 


(LTo C do; £ MD* a (uj,9,E, a)) ’ 

94 where E and a are the electron kinetic energy and local pitch- angle, and D'^ a (w, 8, E, a) 

as is the partial equatorial pitch-angle diffusion coefficient, and the third row shows the 
96 bounce-averaged diffusion coefficients. Note that resonances in come together because 
g? the ui-term can be omitted in the quasilinear resonance condition, u — k\\V\\ — nCl e / 7 = 0, 

98 [e. g., Summers and Thome, 2003], and because the wave spectra are symmetric over 

99 9 = 90°. The “Gauss” lines in Figure 1 show the results for a Gaussian distribution 

100 over x, and reproduce well the equatorial and bounce-averaged diffusion coefficients by 
Albert [2003, Figure 6]. 


Figure 1 


101 
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10 2 Let us first analyze the equatorial pitch-angle diffusion coefficients. For all energies, 

103 Case A is only slightly less than “Gauss” if only |n| — 1 resonances operate, but in the 

10 4 region of |n| > 1 it is about 5 times greater (Figure 1(c) and 1(d), the first row). These 
os dependencies are in good agreement with the previous results by Albert [2003, Figure 10 , 
io« the second row]. For both “Gauss” and Case A, as follows from the second row in the 
10 ? Figure 1 , the contributions from n < 0 are negligible compared to n > 0, especially 
loo for lower electron energies (see Figure 1(a) and 1(b), the second row). Cases B and C 
109 further increase the EMIC wave normal angle, and as a result, suppress the |n| = 1 

no resonances, and for low energies substantially shrink the region of pitch-angles subject 
in to diffusion (see Figure 1(a) and 1(b), the first row). At the same time, they increase by 
ii 2 orders of magnitude the contribution from \n\ > 1 , which operate for greater electron 
n 3 energies, and cover a greater pitch-angle region (see Figure 1(c) and 1(d), the first row). 
H 4 The growing contribution of the n < 0 resonances is more pronounced in Cases B and 
ns C (in comparison to Case A) because EMIC waves become more elliptically polarized 
ns with increasing wave normal angle (see Figure 1, the second row), 
in Overall, in comparison with the field-aligned waves, the intermediate and highly 
lie oblique wave distributions decrease the pitch-angle range subjected to diffusion, and 
ns reduce the equatorial scattering rate by orders of magnitude for low energy electrons 

120 (E < 2 MeV) when only principle |n| = 1 resonances operate. For greater electron 

121 energies, the |n| = 1 resonances operate only in a narrow region at large pitch-angles, 

122 and despite their greater contribution from field-aligned waves, cannot support the local 

123 electron diffusion into the loss cone. In this case, oblique waves operate on the |n| > 1 
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resonances more effectively, and extend the range of pitch-angle diffusion down to the 
loss cone. Note that despite our inclusion of the 7fe + -mode, the above results are in 
qualitative agreement with the results of Glauert and Home [2005, Figures 6 and 7] 
obtained for the equatorial pitch-angle scattering by the 77 + -mode EMIC waves. 

Now we consider the effect of bounce averaging on pitch-angle diffusion coefficients. 
To calculate the bounce-averaged diffusion coefficients, we utilize all the plasma/wave 
parameters used in the above calculation of the equatorial coefficients, and in addition, 
a dipole magnetic field model, and the meridional density distribution from [Khazanov 
et al, 2006]. We further assume that the EMIC waves are confined to mirror points, 
and the wave spectra are equatorial. 

In all considered cases (2), the bounce averaging does not change the shape of the 
diffusion coefficients for energies below 2 MeV (compare the first and third rows in 
Figure 1) but simply reduces the pitch-angle diffusion rates by an order of magnitude. 
For energies 5 and 10 MeV, the peak values of the bounce-averaged diffusion coefficients 
are lower by about a factor of 3 than in the first row of Figure 1. However, the 
bounce-averaged results for E > 2 MeV differ qualitatively from the local coefficients 
for all wave normal distributions. Due to significant scattering at higher latitudes, the 
bounce-averaged diffusion coefficients extend further into the loss cone compared to 
equatorial results. The bounce-averaged results in Figure 1 demonstrate clearly the 
effect of EMIC wave normal angle distribution on relativistic electron scattering. 

Recently, Shprits et al. [2006] showed that the electron lifetime is most sensitive to 
the value of the bounce-averaged scattering rate near the edge of the equatorial loss 
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146 cone, whose value is used to estimate the electron loss timescale [e. g., Summers et al., 

147 2007], Considering the third row in Figures 1(a), 1(b), we can see that the intermediate 
ms and highly oblique wave distributions reduce the scattering rate near the loss cone by 
149 orders of magnitude because only principal |re| — 1 resonances operate. For higher 

iso electron energies (Figures 1(c), 1(d)) when |n| > 1 resonances start to operate, the 
lsi pitch-angle scattering near the edge of the equatorial loss cone depends only slightly 

152 on the wave normal angle distribution, resulting in nearly the same bounce-averaged 

153 diffusion rate for all cases. In other words, there is an electron energy, depending on 

154 specified plasma and/or wave parameters, which separates lower and higher energy 

155 regions with different EMIC wave scattering properties. In the lower energy region, 

156 using a field-aligned wave normal angle distribution leads to a significant overestimate of 

157 the diffusion rate compared to oblique waves. In the higher energy region, the scattering 

158 rate near the edge- of the loss cone almost does not depend on the wave normal angle 
is? distribution. 

iso 3. Local Pitch— Angle Diffusion Coefficient: CRR.ES Based 

161 Calculations 

1 62 3.1. Minimum Resonant Energy 

163 Recently, Meredith et al. [2003] presented an extensive statistical analysis of over 

164 800 EMIC events observed on CRRES to establish whether electron scattering can occur 

165 at geophysically interesting energies (< 2 MeV). In the absence of specific information 
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166 on the wave normal angle, the dispersion relation for strictly field-aligned propagating 

167 EMIC waves was used to obtain the electron resonant energy. For consistency, Meredith 
i6s et al. [2003] included only waves with a high ellipticity (|e| > 0.3) in the survey. 

w9 This yielded a subset of 416 events, the majority of which were identified as L-mode. 
170 Considering only the central wave frequency, u> m , in each wave packet, Meredith et al. 
in [2003] found that in about 11 % of the observations, the electron minimum resonant 

172 energy fell below 2 MeV. These cases were restricted to regions where ui pe /Q. e > 10, and 
to were associated with wave frequencies just below the helium or proton gyrofrequencies. 
174 More recently, trying to increase the above percentage, Loto ’aniu et al. [2006] considered 

173 the entire frequency range for each of 25 EMIC wave packets observed on CRRES during 
ns the initial phase of a geomagnetic storm on 11 August, 1991. These authors also used 
177 the dispersion relation for strictly parallel propagating EMIC waves, and found that in 
17 s comparison with results utilizing LO m only, there are 3 to 4 times more wave packets that 
ns are able to interact with electrons below 2 MeV. 

The minimum resonant energy depends on the wave normal angle, and the 
dependency is stronger in vicinity of the resonant frequencies where the wave number 
grows especially fast. Omitting the w-term in a quasilinear resonance condition 
(cj — &||U|| — nQ e /-f — 0) and taking n = 1, we can obtain the minimum kinetic energy 
required by electrons for cyclotron resonance interaction with EMIC waves, 



iso where E min is the minimum kinetic energy, m e is the electron rest mass, c is the speed 
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iai of light, and k and v are the wave number and electron velocity. Note that equation 

182 ( 5 ) can be obtained from equation ( 7 ) of Summers and Thorne [ 2003 ] by omitting the 

1 83 two smallest terms in their equation. To calculate the electron minimum energy, we 

1 84 select the plasma parameters reported by Loto’aniu et al. [ 2006 , Wave packet # 16 ], 
i 8 s and the results of our calculation are presented in Figure 2 . For 0 = 0 °. as reported 
las in many previous studies [e. g., Summers and Thome, 2003 ], in order to get lower 

187 E m i n , the required wave frequency has to be closer to the He + gyrofrequency (in other 
las words, the wave number should be greater). For most wave normal angles, increasing 
Isa the angle slightly also increases the minimum energy but there is a dramatic decrease 

190 of E min in the region near 8 = 90 °. This transition boundary depends on the wave 

191 frequency. Indeed, there is a resonant wave normal angle (the angle at which the wave 
m number becomes i nfin ite in the “cold plasma” approximation) for any frequency in 

193 the range between and the corresponding bi-ion frequency, and this angle is 

194 closer to 0 = 0 ° if the wave frequency is closer to Q He + . Because of the wave number 

195 increase, the resonant energy decreases dramatically in the vicinity of the resonant 

1 96 wave normal angle, an effect clearly observed in Figure 2 . So in cold plasma, E min is 

197 lower for oblique or highly oblique wave propagation, depending on wave frequency, 
las than for strictly field-aligned propagating EMIC waves. But, of course, the diffusion 

199 coefficient for those wave normal angles should be significant in order to determine the 

200 “physically meaningful” E min , and moreover the cyclotron damping in vicinity of the 

201 He + gyrofrequency can be very strong (see below). 


Figure 2 
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3.2. Pitch-Angle Diffusion Coefficient 

It was demonstrated in section 2 that oblique wave propagation can strongly 
change the effectiveness of both the local and bounce-averaged relativistic electron 
scatterings. At the same time, those results were obtained for plasma parameters 
and wave spectra which were specified independently. So it is important to consider 
the effect of using concurrent observational data. In contrast to section 2, we now 
calculate the local pitch-angle diffusion coefficients using the data for plasma and wave 
parameters reported by Loto’aniu et al [2006]. 

A long duration wave event was observed by CRRES on 11 August, 1991 in the 
interval ~ 0500 — 0700 UT (14.4 - 15.8 MLT) over a magnetic latitude range of —26° to 
—24° and L=6.3-7.6. CRRES was close to apogee in the plasmatrough, and the electron 
density varied slowly from 12 to 17 cm -3 . A total of 25 EMIC wave packets were 
identified both below and above the local He + gyrofrequency \Loto’aniu et al., 2006]. 
In order to estimate the spectral properties of the wave packets, these authors fitted 
a Gaussian distribution to the static wave packet transverse power spectral density. 
Typical FFT data windows and frequency resolutions for the static spectrograms 
were 100 s and 0.02 Hz, respectively. The Gaussian function fit provided the central 
frequencies, u) m , and the spectral semibandwidths, 5u>. The total wave magnetic power, 
5B 2 , was estimated for each wave packet by summing the power spectral density bins in 
the range u> m ± Sea and then multiplying the result by 8lo. Using the full wave spectral 
range, Loto ’aniu et al. [2006] found that electrons with E <2 MeV could interact with 
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223 only three wave packets (16, 17, and 19) if stormtime ion concentration was assumed 

224 (70% H + , 20% He + , and 10% 0 + ). Those packets were the i7e + -mode EMIC waves, 

225 and for the calculation below we selected two of them. The associated plasma and wave 

226 characteristics are summarized in Table 1. Note that to generate this Table we used the 

227 definition of full-width at half maximum (FWHM) as it was given by Loto ’aniu et al 

228 [2006], i. e., FWHM = 2y/2 In 26uj, despite the Gaussian fit ~ exp j— (u — u> m ) 2 

229 Of the packets 16, 17, and 19, wave packet 16 has the most narrow and 19 the widest 

230 distributions, with corresponding power spectral densities presented in Figure 3. 

231 To show the effect of the wave normal angle distribution on relativistic electron 

232 scattering, we use, the wave normal angle distributions ( 2 ), and in addition, a stormtime 

233 ion concentration is assumed. For reference purposes, we also calculate the diffusion 

234 coefficients for strictly parallel/antiparallel propagating EMIC waves. For each wave 

235 packet, the power spectral density is normalized to the corresponding wave magnetic 

236 power 5B 2 shown in Table 1, and this normalization is kept the same for any particular 

237 wave normal angle distribution (2) . In order to estimate the minimum resonant energy 

238 we use yuc from Table 1. For strictly field-aligned wave propagation, as follows from 

239 Figure 2, the energy is about 2 MeV for both wave packets (we can use Figure 2 for 

240 wave packet 19 because cj pe fLl e was nearly the same during both). This minimum 

241 resonant energy exceeds the values presented by Loto’ aniu et al. [2006], especially 

242 for wave packet 16; for this packet and a stormtime ion concentration, they obtained 

243 E min — 0.2 MeV that, as follows from Figure 2, corresponds to a y uc about 0.2496. 

244 Figure 4 shows the results of our calculation for wave packet 16. For strictly 


Table 1 


Figure 3 


Figure 4 
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parallel wave propagation the minimum resonant energy is only slightly below 2 MeV, 
and the diffusion coefficients for field-aligned and intermediate wave propagation are 
only nonzero in Figures 4(c) and 4(d). Cases A and B demonstrate results similar to 
Figures 1(b) and 1(c). Because yuc is very close to the He + gyrofrequency, the minimum 
resonant energy falls below 1 MeV if the wave normal angle exceeds 88°, so that Case C 
may potentially scatter such low energy electrons with an appreciable rate as shown in 
Figures 4(a) and 4(b). Another feature of highly oblique waves is clearly observed in 
Figures 4(d) where the range of pitch-angle diffusion is substantially extended down to 
the loss cone. While Case C exhibits a quite different behavior compared to Figure 1, 
there is a similarity between the diffusion coefficients in Figures 4(d) and 1(c). 

The diffusion coefficients for wave packet 19 are shown in Figure 5. Both Figure 5(c) 
and 5(d) are quite similar, and demonstrate qualitatively the same behavior as in 
Figures 1(a) and 1(b). As follows from Figures 5(a) and 5(b), Case C practically does 
not scatter low energy electrons, mainly because of a lower yuc for wave packet 19 than 
in Figure 4. 


Figure 5 


3.3. Cyclotron Damping Near He + Gyrofrequency and Its Consequence for 
Electron Scattering 

As follows from Table 1, yuc is very close to the local He 1 gyrofrequency 
{VHe+ — 0.25) for both wave packets. In this frequency region, the He + - mode 
experiences strong cyclotron damping due to interaction with thermal He + [e. g., 
Akhiezer et al., 1975]. To demonstrate this, we assume the He + temperature to be 
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T He + — 1 eV, and present in Figure 6 the wave damping rate for the stormtime ion 
composition and plasma parameters observed during wave packet 16. The frequency 
range shown covers approximately the entire wave packet 16. The damping rate for 
yLc has only narrow peak for 0 > 89°, and this region is excluded from the calculation 
of the diffusion coefficients (see equation (2)). For y m , the region of damping near 
ninety degrees extends slightly below 89°, and in addition, small damping appears 
for a field-aligned wave propagation. The situation becomes dramatically different 
for yuc when the jffe + -mode experiences strong damping in the entire wave normal 
angle region; the energy damping rate is 0.5/ r yne+ ~ 7 sec, which is only four times 
greater than the wave period. In all cases, substantial damping takes place only if 
| y — 0.25| ^ k\\V\\ !He +/Q, Ii + 1 where v\\ tHe + is the field-aligned temperature of He + . 
Moreover, we employ a “cold plasma” approximation in our diffusion coefficient 
software (as was done by Loto’aniu et al. [2006]), some must check the validity of this 
approximation. Particularly, the inequality 


Figure 6 


I y - 0.25| » fc, y e+ = e th (6) 

262 should hold. 

263 Inequality (6) is extremely crucial for the diffusion coefficient calculation because 

264 thermal effects should be considered if inequality (6) is violated, but more importantly, 

266 the He + - mode damps strongly in the region \y ~ 0.25 1 For wave packets 16 and 

266 19, inequality (6) is strongly violated in the vicinity of yuc, and waves cannot exist in 

267 these frequency regions, which for T He + = 1 eV, are the ranges £th — 5 x 10~ 3 — 9 x 10~ 2 
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263 and £th — 3 x 10 ~ 3 — 6 x 10~ 2 , respectively. Using these numbers and Table 1, we 

269 conclude that in order to suppress cyclotron damping completely, the He + temperature 

270 should be decreased at least by 1/80 for wave packet 16, and at least by 1/40 for wave 

271 packet 19. Any reasonable change to the temperature assumed in our calculation cannot 

272 eliminate the effect, and can only influence the frequency range subject to cyclotron 

273 damping. 

274 Our conclusion that E3V1IC waves experience strong cyclotron damping near 

275 the He + gyrofrequency contradicts the results of Loto’aniu et al. [2006] because 

276 these authors estimated all their yuc values from CURES data (after filtering, 

277 FFT, and Gaussian approximations). Unfortunately we do not know all the details 

278 regarding data processing used by Loto’aniu et al. [2006], but we know that the wave 

279 frequency resolution in their data was 0.02 Hz. This uncertainty provides the ranges 
™ (.Vlc, Vue ) = (0.20 - 0.25, 0.22 - 0.27) and (y LC , y uc ) = (0.17 - 0.22, 0.22 - 0.27) for 

28 1 wave packets 16 and 19, respectively, that can reconcile our theoretical result with 

282 data reported by Loto’aniu et al. [2006]. So we do not see any reason inequality (6) is 

283 violated, and it must be taken into account. 


284 Let us now recalculate the diffusion coefficients presented in Figure 4, neglecting 

285 contributions from all the partial diffusion coefficients if | y — 0.25| < e t h (keeping all 

286 parameters the same). Note that all the results presented in Figure 1 are still valid 


287 because inequality (6) holds for all those parameters. The results of the recalculation 


288 are presented in Figure 7, and there is a qualitative difference in comparison to Figure 4. 



289 


Now, for all wave normal angle distributions, low energy electron pitch-angle diffusion 
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290 is not possible, and while the 2 MeV diffusion coefficients are nonzero in Figure 7(c), 

291 they are at least partly inside the equatorial loss cone for L as 7.3. For greater electron 

292 energies, the contribution from the high frequency part of the wave power spectral 

293 density decreases. As a result, Figures 7(d) and 4(d) look similar except that diffusion 

294 vanishes at slightly lower pitch-angles in Figure 7(d) than in Figure 4(d), and the 

295 transition between \n\ — 1 and \n\ = 2 resonances is not continued in Figure 7(d) for 

296 Case A. 

297 The results of our recalculation for wave packet 19 are shown in Figure 8. Similar 

298 to wave packet 16, diffusion is not possible for low energies, and Figures 8(d) and 5(d) 

299 are very similar. 

300 In conclusion, we emphasize that as we demonstrated above, the H e + -mode 

301 does not experience significant cyclotron damping by thermal die if y rS Vm (see 

302 Figure 6). So the observed changes in the diffusion coefficients are due to the frequency 

303 region near yuc- and qualitatively correct diffusion coefficients may be obtained by 

304 only considering the region y y m . This result is consistent with the conclusions of 

305 Meredith et al. [2003] regarding the electron minimum resonant energy which were 

306 obtained by considering only the central wave packet frequencies, and suggests that the 

307 number of EMIC wave packets that are able to interact with electrons below 2 MeV 

308 may significantly decrease compared with the estimate of Loto’aniu et al. [2006]. 


Figure 8 
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309 4. Summary and Conclusions 

310 Precipitation of outer RB electrons due to resonant pitch-angle scattering by 

311 EMIC waves is considered to be one of the most important loss mechanisms. The 

312 effectiveness of relativistic electron scattering depends strongly on the EMIC wave 

313 spectral properties, but unrealistic assumptions regarding the wave angular spread 

314 were made in previous theoretical studies. Namely, only strictly field-aligned or quasi 

315 field-aligned waves were considered [Summers and Thorne, 2003 ; Albert, 2003 ; Loto’aniu 

316 et al. 2006 ]. The effect of oblique EMIC waves on relativistic electron scattering 

317 was recently discussed by Glauert and Home [ 2005 ]. For prescribed plasma and wave 
3ie parameters, considering the FT + -mode EMIC waves, they calculated the local diffusion 
sis coefficients and demonstrated that when a realistic angular spread of propagating waves 

320 is taken into account, electron diffusion at ~ 0.5 MeV is only slightly reduced compared 

321 with the assumption of field-aligned propagation, but at ~ 5 MeV, electron diffusion 

322 at pitch-angles near 90 ° is reduced by a factor of 5 and increased by several orders of 

323 magnitude at pitch-angles 30 ° - 80 °. Thus at energies of a few MeV the assumption of 

324 field-aligned wave propagation breaks down, significantly overestimates the pitch-angle 

325 diffusion coefficient at large pitch-angles, and underestimates the local diffusion rate at 

326 smaller pitch-angles by orders of magnitude. 

327 The purpose of the present study was to consider the impact of oblique EMIC 

328 waves on local relativistic electron scattering using simultaneous observations of plasma 

329 and wave parameters from CRRES, and to estimate the effect of bounce averaging. 
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Analyzing 25 EMIC wave packets, and considering the full wave spectral range, 

Loto ’aniu et al. [2006] found that electrons with E < 2 MeV could interact with 
wave packets 16, 17, and 19 only if a stormtime ion concentration is assumed (70% 
H + , 20% He , and 10% 0 + ). Those packets were ffe + -mode EMIC waves, where we 
have selected wave packets 16 and 19 for our analyzes. Results of our study can be 
summarized as follows: 

1. In comparison with the field-aligned waves, the intermediate and highly oblique 
distributions slightly decrease the pitch-angle range subject to diffusion, and reduce the 
local scattering rate by about an order of magnitude at pitch-angles where the principle 
|n| = 1 resonances operate (see Figures 7 and 8). Oblique waves allow the \n\ > 1 
resonances to operate, extending the range of local pitch-angle diffusion down to the 
loss cone, and increasing the diffusion at lower pitch-angles by orders of magnitude (see 
Figures 7(d)). 

2. The local diffusion coefficients based on concurrent plasma/wave parameters from 
CRRES are qualitatively similar to the results obtained for defined plasma parameters 
with model wave spectra (compare Figures 7 and 8 with the first row in Figure 1). So we 
anticipate that the bounce-averaged diffusion coefficients, if estimated from concurrent 
wave/particle data, will exhibit dependencies similar to those we found for the model 
bounce-averaged calculations (see Figure 1, the third row). Those dependencies are: 

3. For low energy electrons, if only principal \n\ = 1 resonances operate, intermediate 
and highly oblique wave distributions (in contrast to field-aligned waves) reduce the 
equatorial pitch-angle range subject to diffusion, and decrease the bounce-averaged 
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352 scattering rate near the edge of the equatorial loss cone by orders of magnitude. This 

353 low energy threshold depends on specified plasma and/or wave parameters, which is 

354 E rs 2 MeV for parameters used in Figure 1 . 

355 4. For greater electron energies, the |n| = 1 resonances operate only in a narrow 
35s region at large pitch- angles (see Figures 1(c) and 1(d)), but due to significant scattering 
357 at higher latitudes, the bounce-averaged diffusion coefficients for field-aligned waves 
353 extend down to the equatorial loss cone. For these energies, oblique waves operating 
359 at |n| > 1 resonances are more effective and provide nearly the same bounce-averaged 
•360 scattering rate in the vicinity of the loss cone as field-aligned waves do (see Figures 1(c) 

361 and 1(d), the third row). 
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Figure 1. Equatorial and bounce-averaged diffusion coefficients versus equatorial pitch- 
angle for scattering relativistic electrons by the /7e + -mode of EMIC waves. Spectral 
parameters and'ion content are given in the text. L=4, and (w pe /Q e ) 2 = it) 3 , where u) pe 

453 and Q e are the equatorial electron plasma frequency and gyrofrequency (without Lorentz 

factor), respectively. The curve “Gauss” is for a wave normal angle distribution adopted 

by Albert [2003]. The second row shows the average resonant number weighted by the 

partial equatorial diffusion coefficient (see the text for definition). 

Figure 2. Minimum resonant energy versus normal angle of the H e + — mode EMIC waves. 

The plasma density and magnetic field are 17 cm -3 and 171 nT, taken from [ Loto’aniu 

454 

et al, 2006, Wave packet # 16]. The ion composition is 70% H + , 20% He + , and 10% 

0 + , and the normalized wave frequency is defined as y = uj/Q h +. 

Figure 3. Transverse power spectral densities for wave packets 16 and 19 obtained by 

455 Loto’aniu et al. [2006]. The solid and dashed vertical lines restrict the frequency range 

± 5u for packets 16 and 19, respectively. 

Figure 4. Local pitch-angle diffusion coefficients for wave packet 16. Calculations are 
based on a stormtime ion composition, 7]h+ = 0.7, rjn e + = 0.2, and rjo+ — 0.1. “W/P 16” 

456 

shows the results for strictly parallel-antiparallel propagating He + modes, and Cases A, 

B, and C are obtained for the corresponding wave normal angle distribution given by (2). 

457 Figure 5. Same as Figure 4, except for wave packet 19. 
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Figure 6. The He + - mode damping rate due to interaction with thermal He + . The phase 
space distribution function for He~ is Maxwellian with T He + — 1 eV, but thermal effects 
are neglected in the real part of the dispersion relation. All other plasma species are 

458 

described in a “cold plasma” approximation. A stormtime ion composition is assumed, 
and the plasma density and magnetic field are taken from [ Loto’aniu et al., 2006, Wave 
packet #16]. 

Figure 7. Same as Figure 4, except inequality | y — 0.25| > k\\v\^ He + /£I h+ is held during 

459 

the diffusion coefficient calculations. 

Figure 8. Same as Figure 5, except inequality \y - 0.25| > k\\V\iH e + /VIh+ is held during 

460 

the diffusion coefficient calculations. 







Wave normal angle (deg) 


64 Figure 2. 





Figure 4. 







0 20 40 60 80 0 20 40 60 80 


Pitch angle (deg) Pitch an g le (deg) 


Figure 7. 





37 


Table 1. Wave Packet and Local Environment Properties Selected From [ Loto’aniu et 


al, 2006] 


Wave 

y-m — 

5y = 

Vlc — 

yuc — 

8B 2 

Bo 

N e 

Packet 

C*W^iT+ 

5oj/ 0. h + 

Vm ~ 5y 

y m + 5y 

nT 2 

nT 

cm -3 

16 

0.23 

0.01 

0.22 

0.24 

2.21 

170.9 

17 

19 



0.22 

0.02 

0.20 

0.24 



0.84 

160.2 

15 







